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A vertically shaken granular medium hosts a blade rotating around a fixed vertical axis, which 
acts as a mesorheological probe. At high densities, independently from the shaking intensity, the 
blade’s dynamics show strong caging effects, marked by transient sub-diffusion and a maximum in 
the velocity power density spectrum (vpds), at a resonant frequency ~ 10 Hz. Interpreting the data 
through a diffusing harmonic cage model allows us to retrieve the elastic constant of the granular 
medium and its collective diffusion coefficient. For high frequencies /, a tail ~ 1// in the vpds 
reveals non-trivial correlations in the intra-cage micro-dynamics. At very long times (larger than 
10 s), a super-diffusive behavior emerges, ballistic in the most extreme cases. Consistently, the 
distribution of slow velocity inversion times r displays a power-law decay, likely due to persistent 
collective fluctuations of the host medium. 

PACS numbers: 45.70.-n,05.40.-a,47.57.Gc 


A corrmrehensive theory of dense granular media is still 
lacking [l|. Unperturbed granular systems may support 
external forces without flowing, behaving like a solid: 
in this configuration, a granular packing may show an 
elastic response to small stresses. Under gentle tapping, 
the granular medium undergoes very slow rearrange¬ 
ments, which resemble the sluggish response of molecular 
glasses [m. The jamming transition, observed when re¬ 
ducing the vibration intensity, or increasing the density, 
is commonly compared to the glass transition in under¬ 
cooled liquids @ . When the energy of the external vibra¬ 
tion is increased or the density is decreased, the granular 
medium enters a liquid-like phase Q which has not re¬ 
ceived as much amount of attention as the glass/solid 
phases or the much more dilute gas phase Q. Neverthe¬ 
less, learning from molecular fluids, the dynamics of the 
liquid phase is rich of information: in particular, it gives 
important hints about the many timescales developing 
when the glass transition is approached from above Q. 
Recent theoretical and experimental insights in the liquid 
phase have highlighted the differences between molecu¬ 
lar glass-formers and driven granular media, in particular 
the presence of superdiffusive behavior [lol4l3j| . 

As rheometers are usually conceived to apply the exci¬ 
tation at the boundaries, standard rheology in granular 
liquids is typically limited by problems of slip and shear 
localization. Other techniques exist to probe the bulk 
response properties of granular media, e.g. magnetic res¬ 
onance imaging. X-ray tomography and high-speed parti¬ 
cle tracking [l^ , confocal microscopy [l^ or multispeckle 


dynamic light scattering [HI- Microrheology [l^ is also 
used, e.g. to probe transient subdiffusive behavior and 
cage effects 17|. However, it is not evident the ability of 
small intruders to probe the collective behavior at large 
spatial scales in the host granular medium. 
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FIG. 1: A: Schematic of the experimental setup. B and C: 
mean squared angular velocity of the blade in the various 
experiments. 

We propose a technique which could be named pas¬ 
sive mesorheology, see Fig. [T]4, inspired from previous 
works The granular medium made of spheres 

of diameter d = 4 mm is placed in a cylindrical container 
of volume ~ 7300 times that of a sphere. The container 
is vertically shaken with a signal whose spectrum is ap¬ 
proximately flat in a range [fmin, fmax] with fmin = 200 
Hz and fmax = 400 Hz. A blade, our probe with cross 
section ~ 8d x 4d, is suspended into the granular medium 
and rotates around a vertical axis. Its angular velocity 
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FIG. 2: Power density spectra of the blade’s angular velocity 
for the two series. A: the number of beads is fixed {N = 2600) 
and the shaking intensity varies. B: the shaking intensity is 
fixed to r = 39.8 and N changes. The frequencies of vibration 
are marked by the yellow bar. Dashed lines are fits of the 
Eq. ([2|l to the experimental data in regions II and III. Thick- 
dashed lines show the limit behaviors f~^ and f~^- 


uj{t) and the traveled angle of rotation 9{t) = 
are measured with a time-resolution of 2 kHz. The blade, 
interacting with the spheres, performs a motion qualita¬ 
tively similar to an angular Brownian motion. Two fam¬ 
ilies of experiments have been performed: a) a series at 
high density {N = 2600), varying the shaking intensity 
r = Zmaxlg G [19.5,39.8], and b) a series at high shak¬ 
ing intensity (T = 39.8), varying N £ [300, 2600] and the 
packing fraction (j) as indicated in the Figures. Fig. [T)3 
and C report the values of the mean squared angular ve¬ 
locity (w^) of the blade in the different experiments. De¬ 
tails on dimensions of the setup and shaken parameters 
are reported in Supplemental Material [2lj. 

Velocity power density spectrum (vpds): in Fig. [2l we 
present our main results in the form of the power den¬ 
sity spectrum of the velocity signal w(t), which is defined 

a® = 2 ntloT We recognize 4 

frequency ranges, denoted as regions I, II, III and IV. In 
the experiments at fixed (maximum) density, N = 2600, 
displayed in Fig. [2]4, the spectrum conserves the same 
qualitative shape, vertically shifted due to the differences 
in F. The most striking properties are observed in re¬ 
gions II and III: S{f) goes from a plateau in region II 


to a roughly parabolic maximum centered at a frequency 
/* in region III. The value of the resonant frequency /* 
slightly shifts from 10 Hz to 20 Hz as F decreases. To 
avoid interference, we ensured that the shaker vibration 
is in a distant region (200 — 400 Hz); we changed such a 
range (including trials with a single frequency, i.e. a har¬ 
monic vibration), obtaining always the same shape S{f) 
with same values of /*. A mechanical resonance is also 
observed at ^ 70 Hz, due to the non-perfect acoustic in¬ 
sulation of the plate on which the couple encoder/blade 
is mounted. In conclusion, the maximum in f* is a res¬ 
onance experienced by the blade in its motion through 
the granular medium, and we interpret it as a transient 
trapping phenomenon, analogous to caging effects in low- 
temperature/high-density liquids. We will see that such 
an interpretation is well supported by other observations 
and by a theoretical model. In both “extremal” regions I 
and IV, S{f) is a decreasing function of /. In particular, 
in region IV, the high-frequency decay presents a power 
law ^ f~^ with 1 < /3 < 2, close to 1 for the lowest values 
of F. This is the evidence of collective effects, in the fast 
processes inside a cage ( as f /*), occurring without a 
characteristic frequency [221] . The decay in region I is also 
anomalous and denotes the emergence of long character¬ 
istic times, possibly larger than the experiment duration 
trOT = 3600 s. Note that the asymptotic diffusion is 
governed by limt_>oo(A 0 ^(t))/t ^ 27r5'(/ —]• 0 ): therefore 
an increasing value of S{f) as / —;• 0 (i.e. at increasing 
time) indicates a superdiffusive behavior, detailed below. 

When the density is reduced, see Fig. 03, the shape of 
S{f) drastically changes. The slope of the decay in region 
I decreases and eventually vanishes: at low densities a 
plateau spans both regions I and H. The maximum in 
region HI is reduced and disappears for N < 1000, called 
the “gas”-phase. The exponent of the power-law decay 
in region IV increases /3 ^ 2. The whole spectrum at 
the lowest densities is well fit by the Lorentzian S{f) = 
;^/[l-|-(27r///7)^], expected for diffusion in diluted gases 
at temperature T with a collision frequency oc 7 [ 2 ^. 
Here T is the probe’s “kinetic temperature” T = I{oP‘). 

Mean-square-displacement (msd). The several phe¬ 
nomena observed in S{f) are reflected in the diffusion 
properties, see Fig. |31 where the mean-square displace¬ 
ment ([A0(t)]^) after a time t is plotted. The four tem¬ 
poral regions corresponding to the frequency regions dis¬ 
cussed above are marked on the graph. Our “cage-like” 
interpretation of the maximum of S{f) in region HI is 
corroborated, in Fig. [SK, by the dramatic slowing-down 
of ([A0(t)]^) in the same region [ 2 ^, resembling the typ¬ 
ical dynamical slowdown in the diffusion of tracers dis¬ 
persed in viscous liquids. At small times (region IV) 
the usual ballistic behavior appears. More remarkable is 
the behavior at large times. All the experiments with 
N = 2600 present a superdiffusive range ^ in region 
H, with a > I. For the largest values of F, this behavior 
changes to a diffusive behavior ^ t in region I. On the 




























3 



t [s] 


FIG. 3: Mean-squared displacement ([A0(t)]^) after a time t 
for the two series. A: N = 2600 and F varies. B: the shaking 
intensity is fixed to r=39.8 and N changes. Black dashed 
lines are guides for the eye. 


contrary, at lower F, the superdiffusive exponent a > 1 
remains the same of region II at large times. In particu¬ 
lar, for r < 31, we find an almost ballistic superdiffusion 
a ^ 2. The situation is very different when the density 
is reduced (Fig. [33): the long time behavior (region I) 
is always of the normal type ^ t. In the most dilute 
cases {N < 1000) the typical scenario of diffusion in a 
gas-like fluid is fully recovered in the form of a mono¬ 
tonic crossover from the ballistic region IV to the normal 
diffusion of region II and I. Changing the size and shape 
of the blade, see [2l|, does not lead to relevant changes 
in the above scenario. 

The study of vpds and msd are consistent with mea¬ 
surements of the velocity-autocorrelation function (vacf) 
C(t), which is the inverse Fourier transform of S{f). In 
the very dilute cases, it is close to a simple exponential 
decay. The most dense and “cold” experiments, on the 
contrary, reveal a vacf with many features: a fast, though 
non-exponential, decay at small times, followed by a 
back-scattering oscillation through negative values, inter¬ 
preted as the “cage”, and finally a slow decay to zero [1^ . 
The final decay of the vacf could also light up the origin 
of superdiffusion. Unfortunately at large times the vacf is 
exceedingly noisy. A more promising way to probe long 
memory effects is to measure the times of persistency, 
i.e. the times during which the signal remains positively 


correlated. Our operative definition consists in two basic 
steps: 1) filter out high frequency oscillations which are 
not relevant for the behavior at large times of A0(t), by 
taking the running average uJsit) = y /j over 

a large time r > Is; and then 2) compute the statistics 
of the times separating two consecutive zeros of uis{t), 
which we call inversion time t^v, see Fig. |4]4 and |4j3 
for experimental samples of ui(t), ujs{t) and Unv The 
statistics of tinv is a natural measurement of long term 
memory of the signal. The experimental probability den¬ 
sity (pdf) of tinv is shown in Fig.|3p-|3F for a few choices 
of parameters and values of r. We observe that it rapidly 
decays with an exponential cut-off smaller or equal than 
^ 10s in all cases where the msd asymptotically showed 
normal diffusion, signaling a finite memory of the dy¬ 
namics. The cut-off apparently jumps to a much larger 
value in the cases where the msd displays superdiffusion: 
in such cases, the pdf decays to zero as a power ^ —2 
or even slower. This is a fair evidence that long memory 
effects arise together with the observed superdiffusion. 
Long memory may be due to a slow rotating creeping 
motion of the surrounding granular medium, which acts 
as a coherent block and drags the blade. Such a “sec¬ 
ondary” motion has long relaxation times due to the in¬ 
volved large inertia. New experiments are being designed 
with the aim of demonstrating this picture. We mention 
that longer experiments (shown in the SM [2l|) still dis¬ 
play superdiffusion that evolves in normal diffusion after 
many hours. 

Modelling the dynamics of the blade. In the most dense 
and cold experiments, the vpds and the msd display 
an intriguing superposition of phenomena over almost 7 
decades of time-scales. It seems hopeless to reduce such 
a complexity to a model with a simple physical inter¬ 
pretation. It is tempting to disentangle the two main 
phenomena, i.e. caging (occurring at times smaller than 
10 “^s) and superdiffusion (more evident at times larger 
than 10s), by decomposing the dynamics into a slow and 
a fast component iij{t) = ojs{t) + ojf (t), with uis{t) defined 
above. Experimental data demonstrate that the stan¬ 
dard deviation oiujfit) is much larger than that of Wg and 
therefore dominates at short times (large frequencies). At 
long times (small frequencies), on the contrary, the fast 
dynamics averages out and uJs{t) emerges as the leading 
signal. For the fast dynamics, we propose a simple inter¬ 
pretation of the transient caging phenomena through a 
“diffusing harmonic cage” (dhc) model: this is a simpli¬ 
fied version of the Itinerant Oscillator model describing 
translational and rotational diffusion of particles in liq¬ 
uids [3 22 1, and reads 


e{t) = ujf{t) 9oit) = \/2Do£,'{t) (la) 

-fA/(t) = -iujf(t) - K[e{t) - 6»o(t)] + v^2^^(t), (lb) 


where ^(t) and ^'{t) are white normal Gaussian noises. 
The model represents the diffusion of a particle in a har- 
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FIG. 4: A and B: angular velocity uj{i) and filtered signal 
ujs{t) (with a running average over a time r = 2 s), for 
N = 2600 and F = 26.8. C,D,E,F: pdf of the inversion times 
calculated (with r = 1 s and r = 2 s) in different experiments. 
From C to F the following configurations are displayed: very 
dilute, intermediate, dense at high energy, dense at low en¬ 
ergy. Colored dashed lines are guides for the eye. 


monic potential with “stiffness” K and unfixed minimum, 
under the effect of a thermal bath at temperature T and 
relaxation time // 7 . The harmonic potential, represent¬ 
ing the cage due to the confining effect of the dense gran¬ 
ular host fluid, is not fixed but moves, as 9o{t) behaves 
as Brownian motion with diffusivity Dq. Motivation for 
this model is twofold: 1) the main features of the vpds, 
i.e. an elastic resonance (region III) and a plateau re¬ 
vealing loss of memory at larger times (region II); 2) in 
the dilute limit it can be rigorously derived [^ . while 
at intermediate densities a series of studies showed that 
memory effects (coming from correlated collisions) are 
well described by a coupling with an additional degree of 
freedom fluctuating at slower time-scales [l^. The vpds 
of the above model can be calculated and reads 


S{f) = - 


1 + iT{2TTff 

212 


TT 7^(27r/)^ [K — /(27r/)^] 


( 2 ) 


Two limiting cases are recovered: when AT = 0, the 
Ornstein-Uhlenbeck process is obtained, with S{f) tak¬ 
ing the Lorentzian form mentioned before. When K > 0 
and Dq = 0, one has the Klein-Kramers process in a fixed 
harmonic potential, and S{f) —0 for / —>• 0, expressing 
the absence of diffusion at large times: the cage does not 
move and fully confines the particle. Formula ([2) fairly 
fits all experimental spectra in regions II and III, with pa¬ 


rameters given in Tables 2 and 3 of [2l| . Reasonably, the 
“cage stiffness” decreases at increasing shaking intensity. 
It also decreases as the density is reduced, and abruptly 
goes to zero for N ^ 1000. The “cage diffusivity” Dq 
rapidly increases with increasing T and with decreasing 
N. A more detailed study of the transition from the cage 
behavior K ^ 0 to the free behavior AT = 0 is postponed 
to future investigations. 

A more ambitious task is to devise a simple mecha¬ 
nism for ujg (t ), leading to superdiffusion. In driven gran¬ 
ular systems it has been observed below the jamming 
transition 11211 and in a few cases above it, where it was 


imputed to “zero”-modes of the host fluid [10|, to Taylor 
dispersion or to turbulence-like cascade effects |^ . 
We stress that, at variance with standard diffusion, for 
anomalous diffusion there is nothing similar to univer¬ 
sality [1^. A systematic derivation of anomalous diffu¬ 
sion is a hard task and it is possible only in few specific 
cases [s^. Of course, the basic ingredient must be an en¬ 
during memory: a way to achieve it is to replace standard 
time derivatives with fractional derivatives in the Fokker- 
Planck equations, an approach which is subject of a vast 
literature 3l|. Such an approach is capable, in princi¬ 
ple, of describing both caging and superdiffusion within 
a single model equation, at the price of losing immediate 
interpretation and plain calculations. In a complemen¬ 
tary approach a coarse-grained value of uig (t) (where 
only the sign of this quantity is traced) follows a contin¬ 
uous time random walk (ctrw): It takes discrete values 
with random transition times extracted from a given dis¬ 
tribution. A simplified version of the ctrw is discussed 
in the SM [2l|. We highlight that the ctrw gets a sug¬ 
gestive experimental confirmation in the observations of 
Fig. |4l indeed the ctrw model quantitatively connects 
the observed slow decay of P{tinv) to superdiffusion. 

Conclusion: an experimental study of dense granular 
mesorheology allows to probe time-scales in a range of 
six orders of magnitude. Such an investigation reveals a 
complex scenario with different dynamical behaviors in 
four frequency or time regions. Three crucial features 
are observed at large density and low temperature: a 
“resonant” caging phenomenon at intermediate scales, 
non-white noise at fast scales, and super-diffusion at long 
times. The caging phenomenon is compatible with a dif¬ 
fusing harmonic cage model, while superdiffusion seems 
to be rooted in the long inversion times appearing in the 
dynamics ujs(t), possibly related to creeping rotating mo¬ 
tion of the granular media. A more detailed investigation 
of transitions is in order: a gas-liquid transition could be 
put in evidence by studying how AT —^ 0 when N de¬ 
creases; a liquid-glass transition could be taking place at 
the lowest values of F or, eventually increasing N > 2600. 
Superdiffusion and large values of Unv could be signaling 
such a transition. Another promising line of investigation 
is active mesorheology, i.e. probing the response of the 
system to the application of an external force, achieved 
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by coupling a controllable motor to the blade (l9l |. 
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SM: SUPPLEMENTAL MATERIAL 


DETAILS OF THE EXPERIMENTAL SETUP 


The scheme of our experimental setup is shown in Fig. lA. A Poly(methyl methacrylate) (PMMA) container 
(diameter 90 mm, maximum height 47 mm and total volume 245 cm^) with an inverse-conical-shaped base, is filled 
with N steel spheres (diameter 4 mm) and mounted on an electrodynamic shaker (LDS V450) which is fed, through 
an amplifier, with a noisy signal whose spectrum is approximately flat in a range [fmin, fmax] with fmin = 200 Hz 
and fmax = 400 Hz. The shaker responds roughly linearly in the same range of frequencies, as verified with an 
accelerometer placed on the container. The elevation of the container at time t is denoted by z{t). We give the 
range of values for the relevant variables: for the maximum acceleration Zmax ~ 20 — 40 g, the maximum velocity 
Zmax ~ 80 — 300 mm/s and the maximum displacement Zmax ^ 0.03 — 0.25 mm. Details are reported below in Table 1. 
A PMMA blade (dimensions 35 x 6 x 15 mm, momentum of inertia / = 353 g-mm^), suspended to a mounting which 
is mechanically isolated from the container/shaker apparatus, is immersed into the granular medium and can rotate 
around a vertical axis placed at the center of the container. The angular velocity uj{t) of the blade and its absolute 
(i.e. without the modulus 27r operation) angle of rotation 6{t) = are measured by an angular encoder 

(AEDA-3300) at a time-resolution of 2 kHz and angular resolution 40000 divisions per revolution. The granular 
medium is in a liquid state because of the vibration. The blade, under the effect of the interactions with the spheres, 
performs a motion qualitatively similar to an angular Brownian motion. It is important to remark that the blade 
can only rotate and therefore moves along the “free flow” direction, i.e. it only probes the fluctuations of granular 
modes parallel to the boundaries of the container. Moreover, the blade’s size is such that, in the dense cases, its 
rotation is only allowed by the rearrangement of a large number of surrounding particles. We have performed two 
families of experiments: experiments at high density {N = 2600) and different normalized (i.e. in g units) shaking 
intensities P = Zmax/g G [19.5,39.8], and experiments at high shaking intensity (P = 39.8) and different densities 
N G [300, 2600]. In Fig. IB and C we have reported the values of the mean squared angular velocity (w^) of the blade 
in the different experiments. 


r 

39.8 

35.2 

30.6 

26.8 

23 

19.5 

Zmin ^max 

0.25-0.06 

0.22-0.05 

0.19-0.05 

0.17-0.04 

0.14-0.04 

0.12-0.03 

Zmin Zmax 

310.4-155.2 

274.5-137.3 

238.6-119.3 

209-105 

179.4-89.7 

152-76 


TABLE I: Equivalence between the normalized shaking intensity T and the physical amplitudes and velocities of the container, 
expressed in mm and mrn.s”^, with fmin ~ 200 Hz and fmax = 400 Hz 


FITS OF THE EXPERIMENTAL DATA: PARAMETERS 


r 

g 

10" 

Do 
-5 s 

,-l 

K/[I{27v 

)^] 

7// 

s~^ 

Tjl 

Tjl 

r/7 

s~^ 

39.8 

46 

± 

4 

128 

± 

3 

121 

± 

3 

0.51 

± 

0.01 



35.2 

31 

± 

3 

151 

± 

4 

134 

± 

4 

0.47 

± 

0.01 

0.53 


30.6 

13 

± 

2 

181 

± 

4 

151 

± 

4 

0.371 

± 

0.009 

0.51 

0.0024 

26.8 

8 

± 

1 

215 

± 

5 

171 

± 

5 

0.307 

± 

0.007 

0.5 

0.0018 

23 

4.6 

± 

0.8 

260 

± 

7 

217 

± 

7 

0.241 

± 

0.006 

0.47 

0.0011 

19.5 

2.6 

± 

0.5 

320 

± 

7 

250 

± 

7 

0.198 

± 

0.005 

0.44 

0.00079 


TABLE H: Fit’s parameters with Eq. 2, series N — 2600 and E € [19.5, 39.8]. 
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N 

Do 

10-3 s-i 

K/[Ii2nf] 

S~^ 

7/7 

S~^ 

Tjl 

s~'^ 

Tjl 

r/7 

S~^ 

2600 

0.46 ± 0.04 

128 ± 3 

121 ± 3 

0.51 ± 0.01 

0.54 

0.0042 

2000 

1.44 ± 0.04 

88 ± 2 

103 ± 2 

0.57 ± 0.02 

0.52 

0.0055 

1300 

7.7 ± 0.1 

84 ± 18 

131 ± 6 

0.97 ± 0.08 

0.54 

0.0074 

700 (gas) 

0 

0 

102 ± 9 

7.0 ± 1. 

1.1 

0.069 

350 (gas) 

0 

0 

79 ± 16 

16.0 ± 2. 

1.1 

0.2 


TABLE III: Fit’s parameters with Eq. 2, series E = 39.8 g and N £ [350, 2600]. 


EXPERIMENTS WITH BLADES OF DIFFERENT SIZES 


We have performed a number of experiments with different blades, changing the total scattering surface and its 
shape. The thickness of the blades being constant, as well as its material, the mass and momentum of inertia are also 
varied accordingly. 

In Figure [SJa we portray the four shapes used. Shape “A” is the one used in the rest of the paper. Panels b and c of 
Figure [5] represent the vpds and the msd respectively, with the four shapes, in the case with F = 30.6 and N = 2600, 
corresponding to packing fraction (j) = 40%. It appears that the general features of the vpds and msd are conserved 
when the size and shape of the blade are changed. Higher or smaller energy is obviously observed, when the inertia is 
reduced (shapes C and D) or increased (shape B). Differences in the mechanical (spurious) resonance at 70 — 100 
Hz are due to some modifications of the mounting mechanism on which the blade is fixed, occurred between the old 
and the new experiments. In the msd of the larger balde (shape B) it is observed - at large times - a saturation of the 
mean squared displacement which we impute to border effects. Note also that at very large time the measurement of 
the msd is noisy and prone to errors. 



a 

A 

35x15 




B 

35x22 


35x5 [ 

D 

18x15 



FIG. 5: Effect of the size and shape of the blade. The four shapes are represented in panel a. In panel b we show the velocity 
power density spectrum for the four cases, with T = 30.6 and N = 2600 which corresponds to packing fraction (j> = 40%. For 
the same parameters, the behavior of the mean squared displacement is displayed in panel c. 


VERY LONG EXPERIMENTS 


The finite size of the experiment suggests that all relaxation times, even if long, are finite: the displacement 
A9, therefore, is expected to reach asymptotically a normal diffusive behavior. We have performed a 12-hours-long 
experiment in order to check such a hypothesis. The mean squared displacement is shown in Fig.[6l Indeed, at large 
times the superdiffusive behavior shows a crossover toward normal diffusion, which is slowly approached after many 
hours. 
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FIG. 6: Mean squared displacement in a long experiment (total time 12 hours), with N = 2600 and F = 26.8. 


THE CONTINUOUS TIME RANDOM WALK MODEL 


A possible mechanism for superdiffusion is the continuous time random walk model (ctrw) [30l . l32 |. In the ctrw 
the variable (here the velocity) performs a random walk taking discrete values LOf. it remains constant for a random 
time t and then jumps to a new value. The choice of the statistics (pdf) of t is the main parameter of the model and 
determines the properties of the autocorrelation of the velocity and, therefore, of the mean squared displacement. 

Here we consider in details a simplified version of the ctrw model, where the velocity can take only two values. Such 
a simplification does not change the substance of the demonstration, which can be easily generalised to models with 
many possible values. Note that a comparison between this model and the experimental results should be made only at 
long time scales (i.e. such that A0(t) is well approximated by ujs(s)ds) and in the spirit of a strong coarse-graining, 
as if a dramatic loss of resolution allows to measure only the sign of the velocity. 

In the model, ojs changes its value at times with n € [0, m] and t^ = t. In the interval between and tn, ujs 
takes the value wo,n € [“Wo,a;o]. The intervals tinv,n = tn — tn-i, with tinv,n = t, are random and independent 
with the following probability density function 


g{x) a;e[0,U] 

P(tinv = x) X ^ X G {t*,tjnax] i 
0 X > tjnax- 


( 3 ) 


where g{x) is a smooth function: its shape and the value of t* (assumed to be always finite) are not relevant for the 
long time behavior of the displacement. The cut-off tmax may eventually be taken to oo (see below). 

For the msd therefore we have 


([A 0 (t)] 2 ) = 



— ^0 ^ '^, {tinv,n) T ^^{x>0,n^0,k) {tinv,ntinv,k) 
n^k 


( 4 ) 


n—1 


and the values being independent and with zero average, one may drop the last term obtaining 

m 

{[Ae{t)]^) = wg = rnu}l{tLv,n) 

1 

The presence of the cut-off allows us to write 


( 5 ) 


\^inv/ — 

m 


( 6 ) 


and therefore one finally has 


imt)?) 


= Wn 


:(0 

' (t- ) 

\^inv / 


t 


( 7 ) 
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The dependence on the cut-off (at large values of tmax) of the moments of order fc > 1 of P{tinv) are determined 
by the distribution at times larger than t*, i.e. 


, lim ^ lim 

tmax—>00 


const. g > fc -I- 1 

9<k+l, 


x'^-Ux 

The last equality, together with Eq. 0 , gives three possible behaviors for the msd at large times t>tr, 


( 8 ) 


• for g > 3: normal diffusion, independently of the cut-off tmax', 

• for 2 < g < 3: ([A6»(t)]2) - t ; 

• for 1 < g < 2 : ([A 0 (t)]^) ^ 

^max 

For 1 < g < 3, in the range t < tmax, one expects ([A0(t)]^) . The value of /3 can be determined by asking 

that it matches, at t tmax, the asymptotic behavior given above: for instance in the case 2 < g < 3 the behaviors 
([A0(t)]^) ^ t^ax ^ ([A^(t)]^) ~ t^ can match only if ,5 = 4 — g; the same argument for 1 < g < 2 gives (3 = 2. 

The situation is therefore summarized here 


([A0(t)]^) ^ 


t g > 3 or t > tmax 

t'i-g 2 < g < 3 and t < tmax 

t^ 3 < g < 2 and t < tmax- 


( 9 ) 


Fig. 4 of the Letter shows that P{tinv) has an exponential cut-off at tmax ^ 10 seconds for all the experiments where 
diffusive behavior is observed, that is in the most dilute or energetic cases. The cut-off rapidly grows and reaches 
times of the order of the duration of the experiment when N = 2600 and T <31: a power law tail with g < 2 emerges. 
Equation ([S]) makes this power law decay fully consistent with the ballistic behavior observed in Fig. 3 of the Letter. 



